arXiv:1503.02439v3 [astro-ph.CO] 31 Aug 2015 


Cosmic constraint on the unified model of dark sectors with or without a cosmic 
string fluid in the varying gravitational constant theory 
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Observations indicate that most of the universal matter are invisible and the gravitational constant 
G(t) maybe depends on the time. A theory of the variational G (VG) is explored in this paper, 
with naturally producing the useful dark components in universe. We utilize the observational 
data: lookback time data, model-independent gamma ray bursts, growth function of matter linear 
perturbations, type la supernovae data with systematic errors, CMB and BAO to restrict the unified 
model (UM) of dark components in VG theory. Using the best-fit values of parameters with the 
covariance matrix, constraints on the variation of G are (^■) z=3 . 5 ~ 1 . 0015 tg 0075 an d (§ )tod ay — 

— 0 . 7252 l 2'|645 x 10 _13 j/r , the small uncertainties around constants. Limit on the equation of 

state of dark matter is wo dm = 0.0072+00170 with assuming wode = — 1 in unified model, and dark 
energy is wode = —0.9986toooii with assuming wo dm = 0 at prior. Restriction on UM parameters 
are B s = 0.7442+00)3^+00292 an d « = 0-00021q;q209-o’o422 with la and 2cr confidence level. In 
addition, the effect of a cosmic string fluid on unified model in VG theory are investigated. In this 
case it is found that the ACDM (f 2 s = 0, /3 = 0 and a = 0) is included in this VG-UM model at la 
confidence level, and the larger errors are given: f l s = -O.OlO6lQ 03J5to.0509 (dimensionless energy 
density of cosmic string), (^) z=3 . 5 ~ 1.00081J;053° and ( %)tod ay ^ -0.349612^3135 x 10 -13 i/r -1 . 


PACS numbers: 98.80.-k 

Keywords: Time-varying gravitational constant; unified model of dark components; equation of state (EoS); 
structure formation. 


I. Introduction 


Gravity theories are usually studied with an assumption that Newton gravity constant G is a constant. But some 
observations hint that G maybe depends on the time lj| , such as observations from white dwarf star Q, Q, pulsar [ij] , 
supernovae jf| and neutron star jfj . In addition, cosmic observations predict that about 95% of the universal matter 
is invisible, including dark matter (DM) and dark energy (DE). The unified models of two unknown dark sectors (DM 
and DE) have been studied in several theories, e.g. in the standard cosmology [rj-Ql, in the Horava-Lifshitz gravity 
10(, in the RS 11] and the KK higher-dimension gravity 12|. In this paper, we study the unified model of dark 
components in theory of varying gravitational constant (VG). The attractive point of this model is that the variation 
of G could result to the invisible components in universe, by relating the Lagrangian quantity of the generalized 
Born-Infeld theory to the VG theory. One source of DM and DE is introduced. In addition, give n that cosmic string 


have been studied in some fields, such as in emergent universe 


13l [l4j, in modified gravity 13 , in inflation theory 
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H, 


and so on 


lTj-^l- Here we discuss the effect of a cosmic string fluid on cosmic parameters in VG theory. Using 


the Markov Chain Monte Carlo (MCMC) method [2l|, the cosmic constraints on unified model of DM and DE with 
(or without) a cosmic string fluid are performed in the framework of time-varying gravitational constant. The used 


cosmic data include the lookback time (LT) data [22), [23], the model-independent gamma ray bursts (GRBs) data 24(, 


the growth function (GF) of matter linear perturbations 


errors 


the type la supernovae (SNIa) data with systematic 
33j, the cosmic microwave b ackg round (CMB) (34j . and the baryon acoustic oscillation (BAO) data including 


the radial BAO scale measurement [351 ] and the peak-positions measurement 


36-38). 


II. A time-varying gravitational constant theory with unified dark sectors and a cosmic string fluid 


We adopt the Lagrangian quantity of system 


L = ^ {w) +16 ' £ “) 


(i) 


with a parameterized time-varying gravitational constant G = Goa(t) - ' 3 . t is the cosmic time, a = (1 + z ) _1 is the 
cosmic scale factor, and 2 denotes the cosmic redshift. g is the determinant of metric, R is the Ricci scalar, and 
C u = Cb + C r + Cd + C s corresponds to the Lagrangian density of universal matter including the visible ingredients: 
baryon Cb and radiation C r and the invisible ingredients: dark sectors Cd and cosmic string (CS) fluid C s . Utilizing 
the variational principle, the gravitational field equation can be derived 


R.,, v - ^Rg^ = 8 nGT^ + G(V, t d u G- L - ^V^G" 1 ) 


( 2 ) 


in which R^ is the Ricci tensor, T^„ is the energy-momentum tensor of universal matter that comprise the pressureless 
baryon (Wb = ^ = 0), the positive-pressure photon ( w r = = |) , the CS fluid ( w s = = — 3 ) and the unknown 


Pr 


dark components (wd = —). w is the equation of state (EoS), p is the pressure and p denotes the energy density, 

Pd ' 


respectively. Taking the covariant divergence for Eq. ([2]) and utilizing the Bianchi identity result to 


or its equivalent form 


3 #(^) 2 + 3-§ + 8n[Gp + Gp + 3HG(p + p)] = 0 

Cr Qj (jr 


( 3 ) 


3iL/3[(/3 - l).ff 2 -H]+ 8 t iG[p + 3 H(p + p)- /3Hp] = 0. 

In the Friedmann-Robertson-Walker geometry, the evolutional equations of universe in VG theory are 

87rGn 


H 2 = 


-a-Pp-PH 2 


( 4 ) 


2—b H 2 = - 8 n G 0 a~Pp - pH 2 - p 2 H 2 - P~. 
a a 


( 5 ) 


From Eq.(U), we can see that a CS fluid can be equivalent to a curvature term in constant-G theory, while this fluid 
could not be equivalent to the curvature term in the VG theory due to the term a -/3 multiplying the density. Combing 
the Eqs. ([31), ([3]) and ([5]), we have 

2 


p + 3 H (p + 


2 + p 


20 P ~ P 2 „ 

P) = ^ , Q Hp. 


2 + p 


( 6 ) 
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’’Dot” represents the derivative with respect to cosmic time t. Integrating Eq. (JG)) can gain the energy density of 

— / 3 2 — 2/3 — 6 — / 3 2 — 4/3 — 8 

baryon pb oc a ^ , the energy density of radiation p r oc a s+P and the energy density of cosmic string 

— / 3 2 — 4 

p s oc a 2 +p . Relative to the constant-G theory, the evolutional equations of energy densities are obviously modified 
in VG theory for the existence of VG parameter /?. 

We concentrate on the Lagrangian density of dark components with the form Cd = [ 1 — (V 

from the generalized Born-Infeld theory jdij], in which V(<p) is the potential. Relating this scalar field <p with the 
time-varying gravitational constant by p(t) = G(f) _1 , it is then found that the dark ingredients can be induced by 
the variation of G. The energy density of dark fluid in VG frame complies with 

Pd = Po d [ B s + (1 - (7) 


here parameter /? reflects the variation of G, a and B s = g ^ 2 B +6 i+T 4 - are m °del parameters. Eq. 0 shows that 

p + / p 0VG _ GCG 

-i 1 Q_|_ P — 

the behavior of pa is like cold DM at early time 1 (for a <C 1, p d ~ pod(l — B s ) T +°a h 2 +t* ), and like cosmological- 

1 

constant type DE at late time (for a> 1, pd « podBs +a ). Then Eq. Q introduces a unified model (UM) of dark 
sectors in VG theory (called VG-UM). The Hubble parameter H in the VG-UM model reads 



H = 


/ Q , 2/3 2 + /3w 1 . n I —2/S 2 —4/3 — 6 -2S 2 -6S-8 -2S 2 -2S-4 

{D od [R s a _ ^( 1+ “) + (1 - B s )a ( 3+ 2 +/< )(!+“)] i+« + Q b a 5+3 + fl r a 5+3 + fl s a 


- 2 / 3 2 - 4 / 3-6 


2 + S 


( 8 ) 


with Hubble constant H 0 and dimensionless energy densities fib = 87T ^^ 0b , fl r = 8n ^par ^ ^ _ s-kG^pos ^ anc [ 
+ fib + + fl s = 1 + /3. For /3 = 0, bove equations are reduced to the standard forms in the constant-G theory. 


III. Data fitting 


A. Lookback Time 


Refs. 


41 


42( | define the LT as the difference between the current age to of universe at z = 0 and the age t z of a 


light ray emitted at z, 


t L (z) = f 
J 0 


dz 


(9) 


/o (1 + z')H(z'Y 

Then the age t{z{) of an object at redshift z t can be expressed by the difference between the age of universe at Zi and 
the age of universe at zf (object was born) [221 ]. 

. , [°° dz r°° dz 

t(Zi ) = J 


(1 + z')H(z') J ZF (1 + z')H(z) 
For an object at redshift Zi, the observed LT subjects to 


7 T = t L {z F ) - t L (zi)- 


( 10 ) 


tr = tL(z F ) - t L {Zi) = K s - t(Zi)] - K S - t L ( ZF )} = t° bs - t(Zi) - df. 


( 11 ) 


1 


P describes the effect on energy density of dark matter from variation of G. 
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One defines 


Xa 


= E 


[tL^-t^^df)} 2 , [to -Cl 2 


( 12 ) 


with af obs + of = erf,. a t cbs is the uncertainty of the total universal age, and <Ji is the uncertainty of the LT of galaxy 
i. Marginalizing the ’nuisance’ parameter df results to 


431 


poo d2 

X 2 lt(Ps) = -2 ln j d ( d f ) exp(-Xa Se /2) = A - -jj + 


[to — t?) S ] 2 rt r / K r / & M 

~ZZ - 21,,| V 2c erfci vW )] ’ 


(13) 


where A = Ar,-B = JU and A = tL(zi) — [tg hs — t(zi)], respectively. p s denotes the theoretical 

model parameters, erfc(ir) = l-erf(a;) is the complementary error function of x. The observational universal age at 


4-obs 


Zi 

0.10 

0.25 

0.60 

0.70 

0.80 

1.27 

0.1171 

0.1174 

0.222 

0.2311 

0.3559 

0.452 

0.575 

0.644 

0.676 

0.833 

0.836 

0.922 

1.179 

ti 

10.65 

8.89 

4.53 

3.93 

3.41 

1.60 

10.2 

10.0 

9.0 

9.0 

7.6 

6.8 

7.0 

6.0 

6.0 

6.0 

5.8 

5.5 

4.6 

Zi 

1.222 

1.224 

1.225 

1.226 

1.34 

1.38 

1.383 

1.396 

1.43 

1.45 

1.488 

1.49 

1.493 

1.51 

1.55 

1.576 

1.642 

1.725 

1.845 

ti 

3.5 

4.3 

3.5 

3.5 

3.4 

3.5 

3.5 

3.6 

3.2 

3.2 

3.0 

3.6 

3.2 

2.8 

3.0 

2.5 

3.0 

2.6 

2.5 


TABLE I: The 38 data points of galaxy age [22|, |23j]. The first 6 data are from Ref. 


■Q. 


B. Gamma Ray Bursts 


In GRBs observation, the famous Amati’s correlation is log ■ 


= a 


b log 


300 keV 


45i 


46J, where E iso = 


Airdf Ssolo /(1 + z) and E p ^ = E Pi0 b s (l + z) are the isotropic energy and the cosmological rest-franre spectral peak 
energy, respectively, is the luminosity distance and Sboio is the bolometric fluence of GRBs. Ref. 
a mo del-independent quantity of distance measurement, 


47| introduced 


- / ^ v p( z ) ( \ (!' 

r p{ z i) = — s,r P (z) = — 


>) 1/2 Ho , , , , d L (z) 


r P (z 0 ) : 


-r(z),r{z) = 


1 + z 


(14) 


with zq being the lowest GRBs redshift. For GRBs constraint, Xgrbs ^ as a f° rm 


Xgrbs(Ps ) = [A r p {zi)\ ■ ( Cov G l RBs )ij ■ [A r p {zi)\ 


(15) 


in which A r p {zi) = r p a (zi) — r p (zi ), and ( Cov GRBs )ij is the covariance matrix. Using 109 GRBs data, Ref. [241 
obtained 5 model-independent datapoints listed in table El where a(r p (zi)) + and cr(r p {zi)) are the ler errors. The 
{fp(zi)} correlation matrix is Q] 


( COVGRB ) 


/ 1.0000 

0.7780 

0.8095 

0.6777 

0.4661 \ 

0.7780 

1.0000 

0.7260 

0.6712 

0.3880 

0.8095 

0.7260 

1.0000 

0.6046 

0.5032 

0.6777 

0.6712 

0.6046 

1.0000 

0.1557 

^ 0.4661 

0.3880 

0.5032 

0.1557 

1.0000 ) 


(16) 


















































5 


with the covariance matrix 


(CovGRB)ij = a{r p {zi))a{r p (zj)){Cov G RB)iji (17) 

where a(r p (zi)) = a(r p (zi)) + , if r p (z) > r p {z ) data ; a(r p {zi )) = a(r p {zi))~ , if r p (z) < r p (z) data . 


Number 

z 

r data ( z ) 

cr ( r p ( zi )) + 

a ( r p ( zi ))- 

0 

0.0331 

1.0000 

— 

— 

1 

1.0000 

0.9320 

0.1711 

0.1720 

2 

2.0700 

0.9180 

0.1720 

0.1718 

3 

3.0000 

0.7795 

0.1630 

0.1629 

4 

4.0480 

0.7652 

0.1936 

0.1939 

5 

8.1000 

1.1475 

0.4297 

0.4389 


TABLE II: Distances calculated by using the 109 GRBs data via Amati’s correlation 



C. Growth Function of Matter Linear Perturbations 


The Xgf can be constructed by the growth function of matter linear perturbations / 


Xgf(Ps) — y 


[fth {Ps 5 + ) fobs(Zi )] 


(18) 


where the used observational values of f 0 b s are listed in table Hill / is defined via /(a) = a , with D = y ■ 

' denotes derivative with respect to a . So in theory, / can be gained by solving the following differential equation in 
VG theory 

tt 2 c 


, M (a) 1 

D (a) + . . + - 

U ( a ) a 


4 + 2/3 + 2/3 2 / 6 + 2 /3 +/3' 2 -e-w-n* 

-\D (a)- /n - ;j^ N2 „2 a 5T3 D y a > = b- 


(19) 


a(2 + /3) J w (2 + /3) 2 H(a) 2 a 2 

For /3 = 0, above equation reduces to the constant-G theory. The derivation of evolutional equation D(a) in VG 
theory are shown in appendix. Comparing with the most popular ACDM model, the effective current matter density 
can be written, flom = fit, + (1 + /3 — fl s — fib — fl r )(l — B s ) for VG-UM. Obviously, for /3 = 0 it is consistent with 


the form of flgm in UM of constant-G theory 


48-5' 


Zi 

0.15 

0.22 

0.32 

0.35 

0.41 

0.55 

0.60 

0.77 

0.78 

1.4 

fobs 

0.51 ± 0.11 

0.60 + 0.10 

0.654 + 0.18 

0.70 + 0.18 

0.50 ± 0.07 

0.75 + 0.18 

0.73 + 0.07 

0.91 + 0.36 

0.70 ± 0.08 

0.90 + 0.24 

Ref. 

[25, 26] 

[27] 

[28] 

[29] 

[27] 

[30] 

[27] 

m 

[27] 

m 


TABLE III: The observational data of growth function fobs- 
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D. Type la Supernovae 


We use the Union2 dataset of SNIa published in Ref. Q- In VG theory, the theoretical distance mod¬ 
ulus Hth(z) is written as fi t h(z) = 5 log 10 [£> L (x)] + ^ log 10 ^ + p 0 , where D L (z ) = ^(1 + z) 2 D A (z) and 

H~ i 1 

fj ,o = 5^°5io(jv^:) + 25 = 42.38 — hlogwh. h is a re-normalized quantity defined by H 0 = 100 h km s Mpc” 

Da{z) = 


=sinn[-y/|f2fe| / 0 ~ ^ z ,, ] is the proper angular diameter distance, here sinn(yj?2fcja;) denotes 


_ (l-KVl^l LV I--VI Jo H{z')\ 

sin(1SI*,|a;), ^\Vtk\x and sinh(y / |flfc|x) for f Ik < 0, fl*, = 0 and f h > 0, respectively. Cosmic constraint from 


51 


SNIa observation can be done by a calculation on 

2 / {Vth{Ps,Zi)~ Li obs ( Zi )} 2 _ ^ {5log w [D L {p s ,z)} + ^\og w ^-m obs (zi) + M'} 2 

XSNIa\Ps) ~ 2^ ^2 — ^2 ' 

SNIa SNIa i 

where p 0 bs{zi) = m 0 bs(zi) — M is the observed distance moduli, with the absolute magnitude M. The nuisance 
parameter M =_Uo + M can be marginalized over analytically, XsNiaiPs ) = —2 In exp \— ^XsNia(Ps> At')] dM’ 
resulting to 


61- 


X 2 SNIa(Ps) = A-(B 2 /C), 


( 21 ) 


where 


15 G 15 G 

A = E "L 5lo Si 0 lD L (ps,Zi)\ + — log 10 — - m obs (zj)} ■ C~ x ■ {5log 10 [D L (p a ,Zj)] + — log 10 — - m 0 b s (zj)} 

SNIa 


Go 

_^ 25 G 

B = E C~ x ■ {5log 10 [D l ( p s ,Zj)\ + — log 10 — - mn obs {zj)} 

SNIa 4 G ° 


C= E c * x - 


( 22 ) 


SNIa 


The inverse of covariance matrix C A with systematic errors can be found in Refs. 


33, 


7lJ. 


E. Cosmic Microwave Background 


XCMB haS a form 



Xcmb(.Ps ) = A di[Cov 1 (d i {p s ),d j (p s ))][Ad i ] t , 


(23) 


with A di{p s ) = df eory {p s ) - d° bs . 9-year WMAP gives d° bs 
the corresponding inverse covariance matrix 


[Ia{z*) = 302.04, R(z*) = 1.7246, z* = 1090.88], and 


Cov 1 


/ 3.182 18.253 - 1.429 \ 

18.253 11887.879 - 193.808 

V -1.429 - 193.808 4.556 y 


(24) 


z* = 1048 [l + 0.00124(f2f,/i 2 ) °- 738 ] [l + gi(flo m h 2 ) 92 ] is the redshift at decoupling epoch of photons with g\ = 
0.0783(f2&/i 2 ) -0 - 238 (l + 39.5(flb/i 2 ) 0 - 763 ) -1 and g 2 = 0.560 (l + 21.l(O fa /i 2 ) 181 ) _1 , l A (p,;z *) = (1 + is 

the acoustic scale, and R{jp s \z *) = + z*)D A (p s ; z*)/c is the CMB shift parameter. 
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F. Baryon Acoustic Oscillation 


The radial (line-of-sight) BAO scale measurement from galaxy power spectra can be depicted by 

a 


( 25 ) 


Two observational values are A zbao{z = 0.24) = 0.0407 ± 0.0011 and A Zbao{z = 0.43) = 0.0442 ± 0.0015, 

c s is the sound speed of the 


35]. Here r s (z) is the comoving sound horizon size r s = c f c 


s dt 


0 a 


respectively 

photon-baryon fluid, c~ 2 = 3 + § x (n^y)a. z d denotes the drag epoch, z d = i+os^ 9 (noJh'^ ^ I 1 + h(^bh 2 ) b2 } 
with b\ = O.313(flo m /i 2 ) _0 - 419 [l + 0.607(flomft 2 ) 0 ' 674 ] and 62 = 0.238(flo m /i 2 ) 0 ' 223 . 

The measurement of BAO peak positions can be performed by the WiggleZ Dark Energy Survey 


Degree Field Galaxy Redshift Survey 


37] and the Sloan Digitial Sky Survey 


the Two 
Introducing Dy(z) = [(1 + 


one can exhibit the observational data from BAO peak positions 


X = 



' s \-°d) 

ZV(0.106) 

- 0.336 1 


’ 4444 

0 

0 

0 

0 

0 


r s (zd) 

D v { 0.2) 

0.1905 


0 

30318 

-17312 

0 

0 

0 


r,(z d ) 

Dv (0.35) 

0.1097 

, V~ x = 

0 

-17312 

87046 

0 

0 

0 


r,(z d ) 

Dv (0.44) 

0.0916 


0 

0 

0 

23857 

-22747 

10586 


r s (z d ) 

D v ( 0.6) 

0.0726 


0 

0 

0 

-22747 

128729 

-59907 

V 

r a (z d ) 

0.0592 J 


l 0 

0 

0 

10586 

-59907 

125536 


(26) 


f. H- 


where V 1 is the inverse covariance matrix shown in Ref. 

The Xbao can constructed 

[A z B ao(z = 0.24) - 0.0407] 2 [A z B ao(z = 0.43) - 0.0442] 2 


Xbao(.Ps ) — 


0.0011 2 


0.0015 2 


+ X t V~ 1 X. 


(27) 


X 1 denotes the transpose of X. 


IV. Cosmic constraints on unified model of dark sectors with (or without) a CS fluid in VG theory 

Multiplying the separate likelihoods Lj oc e~ x */ 2 , one can express the joint analysis of y 2 

2 _ 2 1 2 1 2 1 2 1 2 1 2 /oo\ 

X — Xlt + Xgrbs + Xgf + XsNia +Xcmb + Xbao- (28) 

A. The case with a CS fluid 

In order to obtain the stringent constraint on VG theory, we utilize the cosmic data different from Ref. H to 
calculate the joint likelihood. Concretely, the LT data, the GRBs data, the GF data, the SNIa data with systematic 
error and the BAO data from radial measurement are not used in Ref. 0]. After calculation, the 1-dimension 
distribution and the 2-dimension contours of parameters for the VG-UM model with a CS fluid are illustrated in 
Fig. m From Fig. [T| and table IIVI we can see that the restriction on dimensionless energy density of CS is = 
—O.0106l 0 ; 0305 l 0 050 g in the varying-G theory with containing unified dark sectors. In the constant-G theory, one 























FIG. 1: la and 2a contours of parameters for the VG-UM model with a CS fluid (left) and the ACDM (right) model. 



Mean values with limits (VG-UM) 

Best fit (VG-UM) 

Mean values with limits (ACDM) 

Best fit (ACDM) 

— 

0 _ n ni 

— u.u_LUU_o 0305—0.0509 

0.0006 

O _ n rjrjrjo+0.0024+0.0052 

— u.uuuz-o 0024-0.0048 

-0.0004 

p 

n ni 9«+°- 0394 + 0 - 0756 

U.U1ZO_ 0 . 0385 _ 0 .0718 

0.0005 

0 

0 

B s 

n 1A ^'7+0-0147+0.0269 
u -' 1 -0.0145-0.0299 

0.7520 

— 

— 

a 

n n9i£+°- 0757 + 0 - 1504 

U.UZiO o 0781-0.1466 

0.0004 

0 

0 

h 

U.Oyzz Q 0149—0.0289 

0.6981 

H gqi g+0-0100+0-0197 

u.oyiu_ 0 0101 - 0.0193 

0.6930 

100 Q b h 2 

9 oc:on+ 0 - 0555 + 0 - 1154 

Z.ZOOU o .0557-0.1051 

2.2691 

9 OfxOQ+0-0412+0.0815 

Z.ZUOO o 0420—0.0776 

2.266 

^0 de 

H gqoq+O.0165+0.0347 
U.U^OO_ 0 0161-0.0309 

0.7175 

U. luyO_o.0140-0.0294 

0.7126 


TABLE IV: The mean values with limits and the best-fit values of parameters for VG-UM model with a CS fluid. 


knows that a CS fluid with w s = —1/3 is usually equivalent to a curvature term. But, in the VG theory this equivalence 
is lost due to the term a multiplying the density, as shown in Eq.(U). Comparing the VG theory with the constant- 
G theory, it can be seen that the uncertainty of S2 s in VG theory is larger than some results on Gfc in constant-G 
theory. For example, using the same data to constrain other models we have f Ik = —0.0002/)Q;QQ24io!oo48 (with 
model parameter f2 0 de = 0.7098!q;oi4o1o. 0294) i n ACDM model, f Ik = -O.OOOlioioo 25 io .0050 (with model parameters 
B s = O.7665ioioo9s/m!o205 an d a = 0.0209io!oi89ioio373) i n constant-G UM. Taking the ACDM model as a reference, we 
can see that the influence on the fitting value of flk is small from the added parameter B s and a as seen in constant-G 
UM model, while the influence on the value of f2 S is large by the added VG parameter /3 as indicated in VG-UM model. 
From table Hvl one reads VG parameter f) = -O.Oi 28 io o 385 io.on 8 - Other parameters are B s = 0.7457^'oi^io'o^gg 
and a = 0.0216io'.o78iio'.i466- We then find at la confidence level, the flat ACDM model (f2 s = 0, /3 = 0 and a = 0) 
is included in the VG-UM model with a CS fluid. This result in VG theory is same as the popular point that the 
complicated cosmological model is usually degenerate with the ACDM model. 
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B. The case without a CS fluid 



FIG. 2: 


la and 2a contours of parameters for the VG-UM model without a CS fluid (left) and the ACDM (right) model. 



Mean values with limits (VG-UM) 

Best fit (VG-UM) 

Mean values with limits (ACDM) 

Best fit (ACDM) 

p 

Pi rjr\rj , 7 - l - 0 i 0032-(-0.0062 

U.UUU I o.0033-0.0067 

0.0010 

0 

0 

B s 

0 7AA o - 1 - 0-0137+0.0262 

U. .0132-0.0292 

0.7440 

— 

— 

a 

n nnno+ 0 - 0206 + 0 - 0441 

U.UUUZ o 0209—0.0422 

0.0073 

0 

0 

h 

fl c+0-0098+0.0191 

u,uyuj 0.0096-0.0203 

0.6902 

PI <^Q9c:+0.0094+0.0207 
u.uyzo 0 0104—0.0198 

0.6923 

100 Q b h 2 

0 O£17+0-054+0.116 

Z.ZU < —0.051—0.102 

2.256 

9 9^/| 7+0 ■ 0398+0.0789 

1 o.0394-0.0781 

2.262 

^0 de 

0 7rjQO - l - 0.0148+0.0296 

U. 1 UyO_ 0 . 015 o_0.0309 

0.7095 

r» 71 rj 1 +0.0126+0.0270 

U. 1 J-Ul_o.oi35-0.0282 

0.7106 


TABLE V: The mean values with limits and the best-fit values of model parameters for VG-UM model without a CS fluid. 


For the case without a CS fluid, a stringent constraint on VG parameter is (3 = O.0007i 0 ; 0033 l 0 ' 0067 , where a small 
uncertainty at 2a regions for /3 is given. Still, it is shown that the value of /3 is around zero at ler confidence level 
for both cases: including or not including a CS fluid, and the case containing a CS fluid has a larger error for j3 than 
that not containing a CS fluid. In VG theory, the constraint on UM model parameters are B s = 0.7442+oni32io.o292! 
a = 0 . 0002 to'o 2 O 9 -o.O 422 , h = 0 - 6905 -aoo 96 io.o 2 O 3 and 100 fl b h 2 = 2.267±g;g|f±g;io2- At 1(7 confidence level, the value 
of a = 0 is not excluded, which demonstrates that the ACDM model can not be distinguished from VG-UM model 
by the joint cosmic data. Besides the mean values with limits, the best-fit values of VG-UM model parameters are 
determined and exhibited in table El too. As a reference, the ACDM model is calculated by using the combined 
observational data appeared in section III, and the best-fit values and the mean values with limits on ACDM model 
are laid in tableEl In ACDM model, one receives S2 0 de = 0.71011 q;q 73 51q'o282 that is compatible to the effective result 
of D 0 de in VG-UM model. 
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In order to agglomerate and form structure of universe, one knows that the baryonic (and DM) component must 
have a near zero pressure. Given that wb = ^ = ~^ 2 +/y) ~0, ~ 0 or /3 ~ 1 could be solved. From above constraint 

on parameter /3, one can see that the solution (3 ~ 0 is consistent with our fitting result for both cases: including or 
not including a CS fluid in universe. 


V. Behaviors of G with the confidence level in VG-UM theory with or without a CS fluid 



1.01 
1.005 
O |o 1 
0.995 
0.99 

0 12 3 



z 


1.x 10“ 12 


■o |o 


0 


-1.x 10“ 12 


-2.x 1(T 12 



0 12 3 


z 


FIG. 3: The best-fit evolutions of ^ and ^ with their confidence level in VG-UM model containing (or not containing) a CS 
fluid. 



With CS 

Without CS 

{■§^)z=3.5 

1-0008“ 

1.00151™ 

("G ^today 

—O. 3496 I 26.3135 x lO- 13 ^- 1 

-0.7252l 2 ; 3 ®f 5 x lO^V" 1 

(|))z= 3.5 

-1.800±il|;l|| x t0~ 13 yr~ 1 

-3.792l) 2 ;|it x lO^V" 1 


TABLE VI: The best-fit values of ^ and (S with their confidence level in VG-UM model containing (or not containing) a CS 
fluid. 

In VG-UM theory with or without a CS fluid, the best-fit evolutions of (s) with their confidence level (the shadow 
region) are illustrated in Figure [3] by using the best-fit values of model parameters with their covariance matrix. 
’’Dot” denotes the derivative with respect to t. In the VG-UM model with a CS fluid, limit on the variation of 
G at today is (^)today — — 0.3496l26’3i35 x 10 _13 j/r _1 , and at z = 3.5 we have (^) 2=3 . 5 ~ 0.9917io'Q^ and 
(fj), =3 5 ~ — 1.800ti 3 5; 3 g6 x 10“ 13 j/r _1 . For case without a CS fluid, Fig. [3] shows the prediction that the today’s 
value is (§)today — — 0.7252^2 3645 x 10” 13 yr _1 . This restriction on (^)today is more stringent than other results seen 
in table lVlil Also, using the best-fit value of parameter /3 with error the shapes of = (1 + z)P are exhibited. Taking 
high redshift 2 : = 3.5 as another reference points, we find (^A ) 2=35 ~ 1.0015to'oo75 and (^) z = 3.5 — —0.3792ti;2 3 i4 x 
10 _ 12 t/r _1 in the VG-UM model without a CS fluid. It is important to sternly constrain the value of /3, since the 
monotonicity of ^ = —/3H depends on the symbol of /?. Fig. [3] reveal that the behaviors of G and its derivative are 
around the constant-G theory for both cases: including or not including a CS fluid in universe. 
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Observations 

Limits (yr x ) 

Pulsating white dwarf G117-B15A [2] 

| § |< 4.1 x 10“ 10 

Nonradial pulsations of white dwarfs [3J 

-2.5 x 10 -10 < § < 4 x 10 -11 

Millisecond pulsar PSR J0437-4715 [4] 

| § |< 2.3 x HT 11 

Type-la Supernovae [£] 

| < HT 11 

Neutron star masses [6] 

§ = (-0.6 ±4.2) x 10~ 12 

Helioseismology [75] 

| § |< 1.6 x 10" 12 

Lunar laser ranging experiment [76] 

% = (4 ±9) x 10 -13 

Big Bang Nuclei-synthesis [77] 

-3.0 x 10 -13 < § < 4.0 x 10~ 13 


TABLE VII: Limits on the variation of G. 


VI. Behaviors of EoS with the confidence level in VG-UM theory with or without a CS fluid 


The EoS of UM in VG theory is demonstrated 

Pvg-um /3 — 3 


wvg-um(z) = 


B„ 


(29) 


PVG-UM 3 B s + (1-B s )(1 + 2 )( 1 +“)( 3 -«' 

From Figure [I] (left), we can see that Wvg-um ~ 0 (DM) at early time and wvg-um ~ —1 (DE) in the future for 
the VG-UM model with or without a CS fluid. If the dark sectors are thought to be separable, it is interested to 
investigate the properties of both dark components in VG-UM model. Supposing that the behavior of dark matter is 

— / 3 2 — 2/3 — 6 

known i.e. its EoS Wd m = 0 (pdm = Pod m a 2+3 ), the EoS of dark energy in VG-UM model subjects to 


w de = 


Pde 

Pde 


Pvg-um 
Pvg-um — Pdn 


-A 


-X-\-OL _ n n Ot 

PVG-UM PdmPvG-UM 


(30) 


Using the best-fit values of model parameters and the covariance matrix, the evolutions of Wde with confidence level 
in VG-UM model containing (or not containing) a CS fluid are plotted in Fig. [4] (middle). If one deems the behavior 
of dark energy is the cosmological constant i.e. wa = — 1 (pa = — Pa), the EoS of dark matter in VG-UM model obeys 


'UJdm — 


Pdm 

Pdm 


Pvg-um — Pa 
Pvg-um — Pa 


PaPvg-um A 

PVG-UM ~ P^PVG-UM 


(31) 


which is drawn in 


Fig. 


H (right) with the confidence level for two cases (with or without a CS fluid). 



WO VG-UM 

Wodm (with Wode — ~ 1) 

Wode (with Wodm = 0) 

With CS 

-0.7519t°;°i(2 

O.OOOSVq 0304 

—0.9998jfooi25 

Without CS 

-0.7438t°;°(3t 

0.0072t°;°wo 

-0.9986ig;SSu 


TABLE VIII: The best-fit values of wovg-um , wo dm and Wode with their confidence level hinted by VG-UM model with or 
without a CS fluid. 

From Fig. [H we get the current values wodm = 0.0009io!o304 i n VG-UM model with a CS fluid and wodm = 
0.0072+o.oi70 i n VG-UM model without a CS fluid, which have the larger uncertainties than wodm = 0.0010j(o]ooi6 
calculated on the non-unified model of constant-G theory by Ref. ( 7 ^. For the current value u>ode, it approximates 
to -1 with the very small uncertainty for both VG-UM model with a CS fluid (wode = —0.9998f(g'g^25) an d VG-UM 
model without a CS fluid (wode = — 0.9986j(Q'go^). From the best fit evolution in VG-UM model with a CS fluid, we 
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FIG. 4: The evolutions of EoS with confidence level in VG-UM model including (lower) or not including (upper) a CS fluid. 
Evolution of wvg-um(z) (left), evolution of Wde(z) in VG-UM model with assuming Wdm = 0 at prior (middle), and evolution 
of Wdm(z) in VG-UM model with assuming Wde = —1 at prior (right). 


can see that both Wd e {~ —1) and Wd m {~ 0) tends to be constant, but the uncertainties of them are much larger than 
that in model without a CS fluid. For the best-fit evolution in VG-UM model without a CS fluid, Wd e and Wdm are 
variable with the time and Wdm tends to have small deviation from zero (small-positive pressure) at the recent time. 
In addition, at high redshift the uncertainty of Wd e (or Wdm) is enlarged (or narrowed) for both VG-UM model with 
a CS fluid and VG-UM model without a CS fluid. 


VII. Perturbational behaviors in structure formation for VG-UM theory 


The study on the structure formation is necessary for a cosmological theory. We investigate the evolutions of growth 
function / and growth factor D in VG-UM theory. The derivation of evolutionary equation / and D are shown in 
appendix. Using the definition f(a) = S a , we obtain the dynamically evolutionary equation of / 


(l + z )f'~f 2 + (l + z )f- 


E' 4 + 2/3 + 2/3 2 6 -f- 2/3 + /3 2 flo m 


/ + 


6 + 2f3 + /3 2 


E 2 +(3 J ' (2 + /3) 2 E 2 

where prime denotes the derivative with respect to redshift 2 and E{z) = H(z)/Hq. 


(1 + z) 2 +/5 = 0 


(32) 





FIG. 5: The evolutions of f(z ), f(a) and D(a). The solid lines correspond to the ACDM model, the short-dash lines correspond 
to the VG-UM model with a CS fluid, and the dot lines correspond to the VG-UM model without a CS fluid. 

In Fig. El we use the best-fit values of cosmological parameters in table IIVI and IVl to plot the evolutions of growth 
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function / and growth factor D for VG-UM model and ACDM model by numerically solving Eq. m and l|32|) with 
the initial conditions as* = 0.0001, D{at) = ch, D (a*) = 0 and f(cn) = 1. We can see that the evolutions of /(a) for 
VG-UM model (including or not including a CS fluid) fit well as ACDM model, and the behavior of /(*) are well 
consistent with the observational growth data listed in table Hill In VG-UM model with or without a CS fluid, D(a) 
evolves slower (more slow growth of perturbations) than that in the ACDM model. The current value of D(a = 1) in 
ACDM model is approximately 12% larger than that in VG-UM model without a CS fluid. 


VIII. Conclusions 


Observations anticipate that G may be variable and most universal energy density are invisible. The attractive 
properties of this study is that the variation of G naturally results to the invisible components in universe. The VG 
could provide a solution to the originated problem of DM and DE. We apply recently observed data to constrain the 
unified model of dark sectors with or without a CS fluid in the framework of VG theory. Using the LT, the GRBs, the 
GF, the SNIa with systematic error, the CMB from 9-year WMAP and the BAO data from measurement of radial 
and peak positions, uncertainties of VG-UM parameter space are obtained. 

For the case without a cosmic string fluid, constraint on mean value of VG parameter is (3 = O.OOCUtg'gQgglg'Qgg^ 
with a small uncertainty around zero, and restrictions on UM model parameters are B s = 0.7442^0 ^ 32^0 0292 and 
a = 0.0002^Q Q 2 O 9 io 0422 with 1<t and 2cr confidence level. For the case with a cosmic string fluid, restriction on 
dimensionless density parameter of CS fluid is D s = — O-OlOGlg^Uglg glgg in the VG-UM theory. Obviously, the 
uncertainty of is larger than some results on Dj, in the framework of G-constant theory. At lcr confidence level the 
flat ACDM model (f l s = 0, /3 = 0 and a = 0) is included in the VG-UM model. 

Using the best-fit values of VG-UM parameters and their covariance matrix, the limits on today’s value are 
(%)today = — 0-7252+2'3645 x 10~ 13 or (^)today — —0.3496^26'gi35 x 10 _13 yr _1 for the universe with or without a CS 
fluid. And corresponding to these two cases, we finds (^) 2=3 5 ~ 1.0015to;oo75 and (4^-) z=3 . 5 ~ l.OOOSigggg® at 
redshift z = 3.5. If one considers that the DM and the DE could be separable in unified model, EoS of DE and DM are 
discussed by combing with the fitting results. It is shown that wodm = O.OC^tgggyg or wodm = 0.0009+g;o 3 Q4 with as¬ 
suming wode = —1 for VG-UM universe containing or not containing a CS fluid, while there are wode = — 0.9986t o 'gg 3 i[ 
or w 0 de = —0.9998^0 0125 assuming w 0 d m = 0 at prior for VG-UM model with or without a CS fluid. 

Acknowledgments The research work is supported by the National Natural Science Foundation of China 
(11205078,11275035,11175077). 


Appendix A: The growth of structures in linear perturbation theory 

In a sub-horizon region with length scale r < H -1 , the density of DE and cold DM are expressed by p s d e and p s dm > 
respectively. We suppose that DE is not perturbed, DM is perturbed in sub-horizon region. So, we have p s d e = Pde 
for the homogeneous DE in whole universe and p s dm = Pdm + Spdm for the perturbed DM, where pde and pdm denote 
the density of DE and DM in background level, respectively. Obviously, the region of Spdm > 0 will cluster and form 
structure. In analogy to the equation in background level, the evolution of matter density inside the perturbed region 
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can be given by the following conservation equation 


Psdm + 3 h( 


6 + 2/3 + /3 2 
6 + 3/3 


~ Psdm 


2 + 2/3 
2 + /3 


Psdm ) — 0 * 


(Al) 


Symbol ’’tilde” denotes the cosmological quantity in perturbed region. In this region, the local expansion is described 
by h = r/r and the acceleration is 


8 nG(t) n _ „ . /3 2 , 

~\Psdm + Pde + 3 Psdm + 3 Pde) /- 


r 3(2 + /3)^ som Fae (2 + /3) 

which is same as Eq.© for background level. One can define the density contrast of DM 

Psdm 


1 + $ dm — 


Pdr) 


with 5dm > 0. Differentiating Eq. 


with respect to t gives 


(A2) 


(A3) 


• (S —i— 2/? (S —i— 

5dm + (- 2 + f3 - ^dm)(h — H ) + ^ ^ (1 + 5 dm){JlW s dm — Hw dm ) = 0 (A4) 

after using Eqs. (ED and ©■ Taking the time derivative in above equation obtains 

a 2 6 + 2B + B 2 ■ ■ 6 + 6/3 • 

5dm ~ d ™ -h (- " „ ^ )(1 + <W)(/i - -ff) + 0 r, (hw sdrri + - Hw sdm - Hw sdm )(l + <5 dm ) = 0, (A5) 

1 + Odm 2 + (3 2 + [3 

where 

o n f 4 + 2/3 + 2/3 2 tu m „ 3 H 2 t p sdm ^ n ... , 

h H — D ^dm^dm 0 Q (h HjH ( _ 'dJsdm ^dm^dm) (A6) 

2 + /3 2 + /3 2 + (3 p c 

is given by substituting Eqs. © and (IA2I) into H = - — H 2 and h = - — h 2 , respectively. In addition, in calculation 
we used p c = 3H 2 /8nG(t) and ft + H ~ 277. Inserting (IA6I) into (IA5I) results 


(3 2 

— 


1 $ dm 


4 + 2/3 + /3 2 ij ^ ( 6 + 2/3 + /3 2 ^ H 2 Vldm /x , ^2 \ 

-Hodm ( O I o ) o 1 O \5dm + ^dm) 


2 + ^ 


2 + /3 ' 2 + £ 


r 4 + 2/3 + 2/3 2 6 + 6/3 . _ 6 + 6/3 . ■ ■ ■ . 

[ 2 _|_ /j — iJ (hw s dm HWdm) + 2 _|_ \" J 'Mdm “ 1 “ hw s dm Hw^m HWdm ) 

, / 6 + 2/3 + /3 2 ^ 3-ff 2 ( Psdm ~ 0 ... M n I A \ _ o 

+ ( 0 | 0 ) 0 | ~ Wsdm ^ dm ' dJdm )] (1 + 5 dm ) — O. 


(A7) 


2 + /3 2 + /3 p c 

Neglecting square terms of 5 m in ED, we receive the evolutional equation of density contrast in spherical overdense 
region 

- | 4 + 2/3 + 2/3 2 ^ f 6 + 2/3 + /3 2 ^ H 2 Qdm x , r 4 + 2/3 + 2/3 2 6 + 6/3 ^ 

Orfm + 0 n 33 Odm ( 0 q ) ri . n 5dm + [ ri . n 0 n 33 yhW sd m 3lWdm) 


2 + /3 


2 + /3 


2 + /3 


+ 0 0 ( hw s dm + hlbdm H'W s dm 33w s dm ) + ( r, , a )o o ( ” ~ 3v s dm 1 + 5dm) — 0- (A$) 

2 + /3 2 + p 2 + /3 p c 


Taking [3 = 0, the above equation reduces to case of constant G given by reference 
factor D(a), we can rewrite Eq. (1A8I) as follows 


79]. Using the definition of growth 


" E (a) 1 4 + 2/3 + 2/I 2 / 6 + 2/3 +/3 2 D 0 dm ~ 6 ~ 2 , /3 a ~ <32 n/ \ 

D (a) + [ 4-1-—-—— \D (a)-—-— „ a 2 +^ D(a) 


E(a) a a(2 + /3) 


(2 + /3) 2 E(a) 2 a 2 


r (4 + 2/3 + 2/3 2 )(6 + 6/3) ^__ ^ , 6 + 6/3 ^/. r . , ^ 

~\hjWsdm 3iWdm) d~ /rj n\ tt ^sdm + 3l Wdm 33w dm ) 


(2 + j3) 2 a 2 H 


(2 + p)aH 


18 + 6/3 + 3/3 2 


(2 + /3) 2 a 2 p , 


(——Wsdm - ftdmWdm)][l + D(a)5 dm {a = 1)] = 0. 


(A9) 
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The linear regime of cosmological perturbations is valid for all sales during the early radiation dominated era and for 
most sales during the matter dominated era. For uidm — w s dm — 0, above equation reduces to 


D (a) + [ 


E\a) 

E(a ) 


1 

+ - + 
a 


4 + 2/3 + 2/3 2 
a(2 + /3) 


}D\a) 


6 + 2/3 + /3 2 n odm -e-y -/3 2 
(2 + /3) 2 E{a) 2 a 2 ° 


£>(a) = 0. 


Transferring the function from D to f in above equation, we get 


/(«) + 


/ 2 ( a ) 


13(a) 


. 2 2 + /3 + / 3 2 
a 1 (2 + /3) 


)]/(«) 


6 + 2/3 + /3 2 flodm - 6 - o yy g 2 
(2 + /3) 2 £(a) 2 a“ 


(A10) 
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